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Abstract. - In the present work we propose a method to determine fluctuation induced forces 
in non equilibrium systems. These forces are the analogue of the well known Casimir forces, 
which were originally introduced in Quantum Field theory and later extended to the area of 
Critical Phenomena. The procedure starts from the observation that many non equilibrium 
systems exhibit long-range correlations and the associated structure factors diverge in the long 
wavelength limit. The introduction of external bodies into such systems in general modifies the 
spectrum of these fiuctuations and leads to the appearance of a net force between these bodies. 
The mechanism is illustrated by means of a simple example: a reaction diffusion equation with 
random noises. 



Introduction. - In 1948 Casimir [1] noticed that, since in the vacuum the quantum 
electromagnetic field fluctuates rather than vanishing as in the classical case, it would cost 
some energy to introduce a macroscopic body into it. Indeed, the presence of a macroscopic 
body limits the allowed wave-vectors and therefore the energy density. Casimir derived a 
simple expression for the resulting force per unit area acting between two large conducting 
neutral plates facing each other at distance x. The result is: F/A = —hcTT^ / {2AQx^) and the 
presence of he comes from the quantal origin of the force. In the last twenty years, Casimir's 
idea has been extended to some equilibrium systems where fluctuations are of thermal origin, 
rather than quantal and thus the force is proportional to fc^T [2]. The unifying feature of 
systems displaying thermal Casimir forces is the presence of long ranged spatial correlations. 
The best known case are binary fluid mixtures at the critical point [3] , where the Casimir force 
between two plates of area A and separation a; is attractive and proportional to —ksTA/x^. 
The inverse power law character of the force is originated by the divergence of the correlation 
length, ^, as the critical point is approached. The Ornstein-Zernike theory [4] predicts long 
range density fluctuations characterized by a structure factor (defined as the Fourier transform 
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of the density-density correlation function) like: 

where ^ oc |T — Td"^/^. Then the structure factor behaves as S{k) oc at the critical point 
and the free energy of the system scales as kBTA/x'^~^ , from which one derives the above 
expression for the force. 

Casimir forces also appear in equilibrium physical systems where long range correlations 
appear because a continuous symmetry is broken [5]. For instance in liquid crystals, with a 
vector order parameter, the structure factor for a smectic-A phase layered in the ^-direction 
is anisotropic: 

When the fe is parallel to the layers, i.e. k^ = 0, the structure factor diverges as 

Such smcctic liquid crystals also exhibit a Casimir force: F oc —kBTA/x^ [6]. Similar long 
ranged correlations appear also in nematic phases of liquid crystals, or magnetic liquid crystals 
immersed in a magnetic field H [7] , in superfluid films [8] , in solutions of ionic salts or micellar 
systems [9] among many others. 

Although there seems to exist a relation between the presence of long range correlations and 
Casimir forces, this relation has not been elucidated when computing the force. In practice, 
Casimir forces arc obtained by differentiating with respect to the typical distance the free 
energy change [2] induced by the confinement of the fiuctuating field. 

The goal of this Letter is to establish a general relation between correlation functions 
(or the associated structure factors) and Casimir forces. Wc also provide a novel method 
to compute the Casimir force, specially useful for non-equilibrium systems where long range 
correlations are ubiquitous [10] , but arguments based on the free energy cannot be invoked. 

In a previous work [11] wc derived the Casimir force between large intruders in a granular 
system, starting from the structure factors. In that case the balance between the energy dis- 
sipation due to collisions and the energy injection, renders the system statistically stationary, 
but out of equilibrium. As a consequence it develops long range correlations, and the structure 
factors for density, velocity and temperature, as well as the cross structure functions, decay 
as [12]. 

There are a variety of other non-equilibrium systems for which the structure factors have 
been obtained. A major class is represented by systems under spatial gradients, such as an 
isothermal uniformly sheared fluid in the 2;-direction (plane Couette flow) [13, 14]. It was 
foimd that the structure factor in the limit fc ^ behaves as S{k) oc — 7o-%t^, where 70 is the 
shear rate. In a similar fashion, a binary fluid in a Rayleigh-Benard cell with a temperature 
gradient VT below the onset of the instability, shows an enhancement of the structure factors 
proportional to (VT)^/fc'*. Here the fluctuations arc measured in the 3;y-plane. orthogonal to 
the gradient of the temperature [15]. In the same Rayleigh-Benard cell, but heated from above, 
Li et al. showed that the concentration fluctuations in a polymer solution behave as (Vc)^/A;*, 
where Vc is the concentration gradient created by the Sorct effect. Pure concentration gradient 
were studied by Spohn [16] using lattice gases, who showed that the structure factors behave 
like (Vc)^/(A;^ -|- k^), with fco = 7r/(2L), being L the system size. 

Another class of non-equilibrium systems exhibiting long range correlations are those lack- 
ing the detailed balance condition. For instance, it has been shown [17, 18] that systems 
described by the Langevin equation with conservative dynamics and non-conservative noise 
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have long range spatial and temporal correlations. Analogously, non-equilibrium concentra- 
tion fluctuations in reaction diffusion systems can be long ranged under certain conditions [19] 
that will be described in the next section. 

The mesoscopic model. - As an example of how a Casimir force arises in a non-equilibrium 
system with long range correlations, we consider the simple case of a reaction-diffusion sys- 
tem in three dimensions, where the fluctuating density n around the homogeneous reference 
density no obeys the equation 

^ = V-(DV</. + ^i)-A<^ + 6, (3) 

where </> = n — no is the fluctuating field, D is the diffusion coefficient, and A > is the 
relaxation rate. The terms and ^2 describe fluctuations in the diffusive flux and in the 
reaction rate. They are assumed to have white noise spectrum 

{^u{r,mik{r',t')) = T,5i,kS{r-r')6{t-t') 
{Ur,t)Ur',t')) = T2S{r-r')S{t-t') (4) 

where Fi and r2 are the noise intensities. This reaction-diffusion equation models, for instance, 
the set of chemical reactions: A — ^ A + B — ^ 2 A, that do not satisfy the detailed balance 
condition [20]. In this case, (p is the density fluctuation of either the species A or B relative 
to its average stationary density and A = nok2 — ki. When A = the reaction shows a critical 
point at density no = k\/k2- We assume that the pressure in the system is a function of the 
local density p{n). 

Equation (3) not only models the chemical reaction described above, but governs many 
other physical systems, like liquid crystals [21] or superfluid films [2]. All these systems possess 
structure factors similar to those mentioned in the Introduction and the results obtained in 
this Letter can be generalized to those systems. 

The solution of Eq. (3) predicts that, after an initial transient, the density (j) is statistically 
homogeneous and stationary so that in Fourier space we have: 

{(l>k) = 0, {<f>k<P<,) = VS{k)6k,-g (5) 

where the symbol (•) represents the average over the two noises ^1 and ^2- The structure 
factor S{k) is given by (see e.g. Chap. (8.3) of [19]): 



Fi/j^ + Ta _ Fi r/2i:> 
2(L>A;2 + A) "20^ WTk^ 



with F = F2 — FiA/D and fco = \/X/D. The corresponding real space correlation reads 

Fi F 6"^°'" 

^W=2^^(^)+2D^ 

The second contribution, stemming from the /s-dcpendcnt term in (6) , represents fluctuations 
with a correlation length that depends on the reaction parameter. A, and, therefore, are 
of macroscopic size. In particular, near the critical point A — > 0, the correlation length 
diverges. If the reaction satisfies the fluctuation-dissipation theorem [19,22] then Fi — 2kBTD 
and F2 = 2A:bTA, where T is the temperature, implying that F vanishes along with the 
macroscopic correlations. On the contrary, in non-equilibrium systems which violate the 
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fluctuation-dissipation theorem, F does not vanish and macroscopic correlations are present. 
The (5-term in Eq. (7), coming from T\/2D, describes the microscopic self-correlation of the 
particles that a mesoscopic model cannot resolve. Therefore, they can be eliminated. This 
corresponds to subtracting the asymptotic value of S for large values of k. From now on, we 
will consider the macroscopic part (or, equivalcntly, the non-equilibrium part) of the structure 
factor S*{k) = S{k) - limg^oo S{q) = T/[2D{k'^ + k^)]. This is equivalent to suppress the 
vectorial noise and keep only a scalar noise ^ with an intensity F. 

We study now the effect of confining the system between two plates, parallel and infinite 
in the y and ^-directions, located at a; = and x = L, and calculate the force between the 
plates surrounded by a fluid described by Eq. (3). The force derives from the pressure, p{n), 
exerted by the particles over the plates. To proceed, we consider the system in a volume 
L^, periodic in all directions. In this volume we place two plates at distance L with 
non flux boundary conditions at them, as natural for a reacting system. The total volume V 
results divided into two regions: Region I in between the plates of volume L x Ly x Lz, and 
Region II outside the plates of volume {L^ — L) x Ly x Lz- The limit Lx,Ly,Lz — > oo will 
eventually be taken. 

In order to perform the analysis in the two regions let us consider a case of a general volume 
V = X X Ly X Lz, where X = L for Region I and X = L^ — L for Region II. The density field 
is expanded, taking into account the non flux boundary conditions on the a;-direction, as 

0(r, t)^V-'J2J211 Mt) cos(fc,x)e*'=='*'e^^-^^ (8) 

/CgU fcy 

where kx = wUx/X, ky = 2iTny/Ly, kz = 2'Knz/Lz, = 0, 1, 2, . . . and Uy^Uz = ■■■■, —1,0, 1, . . 
The noise ^ is expanded in a similar way with 

{ik{t)Ut'))=lKVT5k,q5{t-t') (9) 

where 5k^q = Sk^^q^Sky,-qySk^,-q^ is a modified 3D Kronecker delta. Moreover the factor 
{ikx = 1/2 if /Cx = and 7fc^ = 1 otherwise) appears because of the non-flux boundary 
condition in the a;-direction. Replacing these expansions in (3) it is found that 

(<^fe) = 0, {Cl>k<l>q) = lkMqVS*{k) (10) 

with the same structure factor S*{k) as in the homogeneous case. Finally, the density field 
fluctuations in real space are given by: 

{<t>{r))=Q, {<t>{rf) = l^V-'Y^^^cos{qxkaxf. (11) 

q 

where q = k/ko and the prime in the sum means that the term Qx = has a factor 1/2. 

The sum in Eq. (11) contains an ultraviolet divergence {q oo). Therefore, in order to 
perform the summation a regularization prescription is needed. We introduce a regularizing 
kernel in Eq. (11) of the form 1/(1 -|- e^q^), that equals 1 for e ^ 0, limit that will be taken 
at the end of the calculations. The election of a rational function instead of an exponential 
one is made to keep the calciilations as simple as possible. The technique of a regularizing 
kernel is equivalent to imposing a cutoff in the g- vectors of the order of qc ~ or to the 
zeta-function regularization method [21]. 

Next, we take the limit Ly,Lz oo allowing us to replace the sums on qy and by 
integrals that can be carried out, with the result: 
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Fig. 1 - Density fluctuations as a function of the dimensionless position kox, for e = 0.05 and 
koL = 1/2. The vertical scale is in arbitrary units and the asymptotic value of the density fluctuations 
has been subtracted. The inset shows the density fluctuations jump at the plate. 



Figure 1 shows the density fluctuations (12) in the x-direction when the plates are located 
at .u = and koL = 1/2, e ~ 0.05 and oo. The set of (7a;-vectors entering in the sum of 

Eq. (12) for Regions I and II arc different. In Region I, the allowed vectors are quantized 
as: = irrix/ {k^L), while in Region II they form a continuum in the limit Lr^. oo. This 
difference produces a jump of {(fr^) at the plate, that is shown in the inset. 

Casimir forces. - As shown in Fig. 1, the density fluctuations in a confined system 
depend on the position but also on the system size, and therefore are different in the region 
in between the plates and the region outside them. If the pressure is a function of the local 
density field p{n), these differences in the density fiuctuations create an unbalance in pressure 
between both sides of each plate, and consequently a net force. To calculate; th{^ force, we 
expand the local pressure around the reference density no, and take the statistical average, 
finding that 

1 d^p 



(j>{x)) = p{no) + 2 ^ 



im')- (13) 



no 



The net force acting on the plate located at a; = i is the difference between the pressure 
inside p{x L~) and the pressure outside p{x L^), that come from density fluctuations 
in Regions I and II respectively. Then the net effective force between the plates by unit area 
is twice the pressure difference in one of the plates: 



167r£l(l - e2) 



1 + 



(14) 



The integral in this equation is simply 27r(e — l)/e, and Eq.(1.431,2) of [23], allows us to 
evaluate the sum. The diverging terms 0{e~^) cancel and the result for the force, in the limit 
of a vanishing cutoff e, is surprisingly simple: 

F/A = F.{l-'^i^) (15) 

where / = k^L^ and Fq = rko{d-^p/dn'^)/{8TrD). Let us note that the final expression of the 
Casimir force, is a universal function of the reduced distance, I = koL^- Moreover, there is 
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no dependence on the cutoff length, as the two divergences in the cutoff, one stemming from 
the discrete sum and another from the integral, exactly cancel each other. The regularizing 
kernel, a technique well known in the field of Casimir forces [24], has allowed us to obtain a 
finite result as a difference of two diverging quantities. 

The analysis of the Casimir force, Eq. (15) can be performed in the limits of far plates 
(/ ^ 1) or near ones {I <C 1). In the first case, (Z ^ 1), implies that the distance Lx ^ k^^ 
and therefore the plates are outside the correlation length, k^^ (see Eq. (7)). Then, one 
expects a very fast decay of the Casimir forces. In the opposite limit (/ <C 1), when the plates 
are well inside the correlation length, the force is much stronger. These forces are 

Ff»^/A = Fo^, F"^"7yl = -Fo^. (16) 

These results are presented in Fig. 2 where we plot the exact force Eq. (15) as a function 
of the dimensionless distance I, together with the far and near plate approximations. 




Fig. 2 - Dimensionless force as a function of the dimensionless distance I, together with the asymptotic 
form for / ^ 1 (dashed Une) and I <C 1 (dotted line). 

Summary and Conclusions. We have calculated the fluctuation induced Casimir force 
between two parallel plates produced by a scalar field whose evolution is given by a reaction- 
diffusion equation. This system is out of equilibrium because the noise terms do not satisfy the 
fluctuation-dissipation relation, developing long range correlations. From these correlations 
one can calculate the pressure that the fluctuating field exerts on the two parallel plates. 

The Casimir force obtained in this way depends on a macroscopic parameter ko that defines 
a characteristic length as ^q"^, which is a combination of the diffusion constant and reaction 
rates, similar to the cases described in other physical systems [16,20]. The calculated force 
exhibits a crossover at plate separation of the order of fc(7^, distance at which the correlation 
have decayed (sec Eq. (7)). For short separations, the correlations play an important role 
and the force shows a strong distance dependence, while for k^^ the density correlations 
have decayed and the force vanishes exponentially. 

Finally the force is repulsive if the pressure is such that d^p/dn^ > 0, as it is usually the 
case for molecular fluids far form the liquid-gas transition. However, the sign might change for 
different geometries of the obstacles (for instance plate-sphere, or two concentric spheres [25]). 
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To conclude, we have developed a method to calculate Casimir forces starting from the 
structure factors. This method is valid for both equilibrium and non equilibrium systems, 
where the usual derivation of Casimir forces based on free energy calculations is not applicable. 
The necessary ingredients to obtain a Casimir force are a fluctuating field with correlations 
that extend over a macroscopic range, and a confinement of the fiuctuation spectrum induced 
by the presence of the plates or other obstacles. 
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